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In magnetically-doped quantum dots changing the carrier occupancy, from open to closed shells, 
leads to qualitatively different forms of carrier-mediated magnetic ordering. While it is common to 
study such nanoscale magnets within a mean field approximation, excluding the spin fluctuations 
can mask important phenomena and lead to spurious thermodynamic phase transitions in small 
magnetic systems. By employing coarse-grained, variational, and Monte Carlo methods on singly 
and doubly occupied quantum dots to include spin fluctuations, we evaluate the relevance of the 
mean field description and distinguish different finite-size scaling in nanoscale magnets. 
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I. INTRODUCTION 


Nanoscale magnets are fascinating systems displaying 
phenomena at the boundary between classical and quan¬ 
tum physics. They reveal important implications for 
fundamental phenomena, such as macroscopic quantum 
tunneling,magnetic polaron formation,tunable 
magnetism,and strongly-correlated states,® as well as 
potential applications in information storage and process¬ 
ing, arising from the superparamagnetic limit,mag¬ 
netic hardening induced by nonmagnetic molecules, 
spin-lasers,and implementations of qubits. 

Despite the significant differences between nanoscale 
magnets and their bulk counterparts, a mean-field de¬ 
scription that could be suitable for bulk magnets in the 
thermodynamic limit, remains also widely used in de¬ 
scribing magnetic ordering in nanostructures. Unfortu¬ 
nately, the appealing simplicity of the mean-field approx¬ 
imation (MFA) can often mask important phenomena. 
Neglecting thermodynamic spin fluctuations can lead to 
spurious thermodynamic phase transitions in small mag¬ 
netic systems. Does that imply that the MFA cannot de¬ 
scribe nanomagnets, or that there are situations in which 
MFA could yield valuable and unexplored insights? 

In this work we show that the applicability of the MFA 
varies between different nanomagnets which also display 
different finite-size scaling and lead to distinct thermo¬ 
dynamic limits. We focus on magnetically-doped semi¬ 
conductor quantum dots (QDs) with the localized impu¬ 
rity spins typically, provided by Mn ions.'^’®d 6-42 q’j^ggg 
systems are multi-carrier generalizations of the magnetic 
polaron formation^^^®^ that can be viewed as a cloud of 
localized impurity spins, aligned through exchange inter¬ 
action with a confined carrier spin. The characteristic 
signatures of magnetic polarons is the presence of high- 
temperature tails in the root mean square magnetization, 
rather than an abrupt vanishing of magnetization at the 
Curie temperature, Tc, expected for bulk magnets.®^ 

To qualitatively distinguish magnetic ordering in dif¬ 
ferent nanomagnets, such as epitaxially-grown or col¬ 
loidal QDs,^’^®“^® we introduce a simple description in 




FIG. 1. (Color online) Qualitative behavior of magnetic or¬ 
dering in nanoscale systems. An example of a system that 
(a) does not and (b) does display a phase transition. The free 
energy, Eq. (1), is shown as a function of the order parameter 
at various temperatures applicable to both (a) and (b). (c) 
The temperature evolution of the order parameter for the free 
energy described in (a) and (b). (d) For high temperature, 
there are two qualitatively different finite-size scalings for the 
normalized order parameter which also extrapolate to distinct 
thermodynamic limits (circles). 


which the relevant free energy functionaU^’^^ is reduced 
to a MFA free energy, F, given as a function of the order 
parameter X, and the absolute temperature, T, 

F{X, T) = go{T) + g, {T)X + + ^-^X\ 

( 1 ) 

where the expansion coefficients, go, gi, a, and 54 are 
functions of T. The quadratic and quartic terms in X 
describe the entropy of the magnetic ions. Compared to 
the conventional Ginzburg-Landau form,®® it is surpris¬ 
ing to see in Eq. ( 1 ) the linear term in X in the absence 
of an applied magnetic field, we describe the origin of this 
term in Sec. III. 
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As shown in Fig. 1 (solid lines), the presence of a lin¬ 
ear contribution in the order parameter, gi ^ 0, has a 
striking consequence. For all relevant the free en¬ 
ergy minimum is attained for a nonvanishing order pa¬ 
rameter [Fig. 1(a)], implying the absence of a spurious 
phase transition. The magnetic order remains finite for 
all T [Fig. 1(c)]. In contrast, for gi = 0 (broken lines) 
Figs. 1(b) and (c) reveal a behavior typically associated 
with bulk ferromagnets displaying a phase transition for 
T > 0. Surprisingly, this simple MFA description pro¬ 
vides two qualitatively different finite-size scaling with 
number of spins (magnetic ions) in Fig. 1(d) which we 
later show accurately reflects the behavior of two classes 
of nanoscale magnets by considering a more rigorous ap¬ 
proach including spin fluctuations. 

After this introduction, in Sec. 11 we provide an 
overview of the employed theoretical methods and discuss 
the importance of the correct choice of the order param¬ 
eter. We then focus on the two classes of nanomagnets 
and the simplest magnetic QD embodiment: (i) a single 
occupancy, implying the finite carrier spin configuration 
of an open shell in Sec. Ill, and (ii) a double occupancy, 
corresponding to the vanishing carrier spin configuration 
of a closed shell in Sec. IV. We explain how these two 
classes of nanomagnets are already qualitatively different 
at the MF level corresponding to (i) gi ^ 0 and (ii) gi = 0 
behavior, in Fig. 1 and how they can be viewed as rep¬ 
resenting magnetic polarons'^’®d6,34,35,55-57 magnetic 

bipolarons,®’^®“®° respectively. We conclude our presen¬ 
tation with the implications for the relevance of MFA to 
nanomagnets and discuss outstanding questions. 


II. THEORETICAL OVERVIEW 

Magnetically-doped QDs are a useful model system to 
study magnetic ordering in nanostructures. Even with 
very different growth techniques (top-town or bottom- 
up), such as epitaxially grown QDs or solution-processed 
colloidal QDs, there are striking similarities in the mani¬ 
festations of their nanoscale magnetism, as well as in the 
limitations of their theoretical description. We illustrate 
different implications of magnetic ordering by focusing on 
(II,Mn)VI QDs, depicted in Fig. 2. These systems dis¬ 
play carrier-mediated magnetism, extensively studied in 
bulk dilute magnetic semiconductors. Since Mn^’*' is iso- 
valent with group II ions, carriers must be created inde¬ 
pendently; for example, excitation of electron-hole pairs 
by interband absorption of light [Fig. 2(a)]. By changing 
the intensity of light can thus change the QD occupancy 
to realize both open- and closed-shell QDs. A realiza¬ 
tion of multiple occupancy in QDs is observed in various 
experiments. 

From the band alignment of conduction and valence 
bands in Fig. 2(b), the ordering of Mn spins located in the 
QD region is dominated by the holes, characterized by 
the Mn-hole exchange coupling /?. The electrons have a 
negligible influence on magnetic ordering. They are spa¬ 



FIG. 2. (Color online) (a) A scheme of QDs grown on a 
two-dimensional wetting layer (WL) in which electron-hole 
pairs are created by inter-band absorption of light. Holes are 
subsequently captured in the QD. Type-II conduction/valence 
band (CB/VB) profile®^ of a II-VI QD doped with Mn spins.®^ 
e and Ecv are the confinement and bandgap energies. 

tially removed from Mn and typically have ~ 5-6 times 
smaller exchange coupling with Mn than holes. 

We first recall the MFA as known from bulk systems, 
but due to its simplicity it is often used in nanomag¬ 
nets where its validity is questionable. A conventional 
mean-field theory is illustrated in Fig. 3, where the most 
probable state of the system is obtained by minimiz¬ 
ing the free energy as a function of the order parame¬ 
ter. For nanoscale systems, one needs to be careful with 
the choice of an order parameter. This often overlooked 
consideration, has striking implications, as depicted in 
Fig. 3(a) and (b), which we further explain on the ex¬ 
amples of magnetic polarons (MPs) and magnetic bipo- 
larons (MBPs), discussed in Secs. Ill and IV. A common 
choice of magnetization ^ as the order parameter, shown 
in Fig. 3(a) leads to spurious thermodynamic phase tran¬ 
sition in nanoscale magnets. In contrast, if an “observ¬ 
able” quantity, such as the exchange energy X, is chosen 
as the order parameter, the phase transition can be re¬ 
moved. This is shown in Fig. 3(b), where there is always 
a minimum in the free energy at all T. Based on the 
choice of the order parameter we can then recover either 
the gi = 0 or the gi ^ 0 behavior of the free energy 
depicted in Figs. 1(a) and (b). 




FIG. 3. (Golor online) (a) The free energy as a function of 
magnetization, reveals a minimum at low-T and vanishes 
at high-T. MFA will give a second-order phase transition, (b) 
Corrected free energy of a nanoscale magnet as a function of 
“observable” exchange energy, X, giving a finite value to the 
order for any finite T. 

Some qualitative trends in magnetic QDs can be ob¬ 
tained from the MFA developed for bulk dilute mag¬ 
netic semiconductors. The exchange coupling of a car- 
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rier (hole) and magnetic impurity (Mn) spins, s and S, 
respectively, can be expressed in terms of two effective 
magnetic fields. The resulting self-consistent equations 
for the average spin (along the direction of the applied 
field or spontaneous Mn magnetization) are different de¬ 
pending on whether the carriers have nondegenerate or 
degenerate distribution. later discuss how 
these two bulk cases have similarities with the MFA ap¬ 
plied to MPs and MBPs. 

For nondegenerate carriers, the self-consistent equa¬ 
tions are,^° 

(,.) = SB. (^(&)) . (2) 

(&) = SB. (^(».)) . (3) 

where Og is the unit cell volume, ric and rii are the carrier 
and magnetic impurity densities, respectively; ks is the 
Boltzmann constant, and Bj{x) is the Brillouin function, 

BAx) = coth (^x) - T ^ (4) 

In the high-T limit of a small ratio of the effective mag¬ 
netic and thermal energies, the expansion in Eq. (3), 

Bj{x) K. X + 0{x^) a: < 1, (5) 

yields a vanishing magnetic response at a critical tem¬ 
perature. 


Unlike many studies of magnetic QDs that do not go 
beyond mean field theory,®’®°’^^“’^"‘ we will utilize two 
methods which include spin fluctuations. The first is 
a coarse-grained approach® in which we discretize the 
QD space into a number of cells, with being the 
number of Mn spins belonging to each grid point where 
projection of the total spin onto the z- 
axis of the Mn contained at the grid point. Within 
a given cell the wave function is slowly varying allowing 
one to neglect the spatial dependence of the carrier and 
spin density. The full partition function is obtained by 
summing over all configurations of the normalized mag¬ 
netization in a given cell. 

The second method is to perform Monte Carlo simu¬ 
lations. Unlike the coarse-grained method, Mn can be 
positioned at many sites allowing for spatial variation 
in the carrier spin density. The Monte Carlo simulation 
seeks approximate solutions to the Schrodinger equation, 
for the Hamiltonian H, at a fixed T, for a finite orthonor¬ 
mal basis 1$) at a given Mn spin {S'z} configuration. The 
calculation entails randomly generating a Mn configura¬ 
tion at a given T and producing a matrix representation 
of H{{Sz}) in a finite basis and solving the eigenvalue 
problem. A metropolis algorithm is used to obtain the 
most probable Mn configuration at a fixed 

Our model to study magnetic ordering in open- and 
closed-shell systems is motivated by the Mn-doped QD 
with type-II band alignment in Fig. 2, was shown ex¬ 
perimentally to support robust MPs.^^ We use the to¬ 
tal QD Hamiltonian, H = He + Hex, with typical two- 
dimensional (2D) nonmagnetic (carrier) and magnetic 
(Mn-hole exchange) parts, where 


fesTc = ^/3al,/ninc^y S{S + l)s(s + 1). (6) 

For the degenerate case, the carrier spin is given by 

(sz) = ^ J def{e) [D{e + sl3alni{Sz)) 

- D{e - sPaQ-riiiSz))] ■ (7) 

where D{e) is the density of hole states. In the high-T 
limit, the integrand in Eq. (7) can be expanded, using 
Eq. (3), and Eq. (5), to give the critical temperature 

kBTc = ( 8 ) 

where fi is the chemical potential. Interestingly, the dif¬ 
ference between the linear and quadratic /^-dependence 
of the Tc in MFA for the two bulk dilute magnetic semi¬ 
conductors in Eq. (6) and (8) is also obtained using the 
MFA for MPs and MBPs, respectively. 

At low-T, where spin fluctuations are small, MFA 
can accurately describe the thermodynamics of a finite- 
size system. However, a careful treatment is needed 
for a higher T, where large spin fluctuations could play 
the dominant role in the thermodynamics of magnetic 

QDs.55.7i 


Nh 


I?® 1 

--V| + 

2m* 2 * 


U, 


Nh^ 


( 9 ) 


h is the Planck’s constant there are Nh holes at the posi¬ 
tion Vi and m* is their effective mass. A harmonic x — y 
conhnement of strength oj is much weaker than the con¬ 
finement along the growth (z) direction, implying effec¬ 
tively 2D system. U is the charging energy. The p — d 
exchange interaction between spins of Mn and confined 
holes has the Ising form^^’^®^™ because of the strong 2 ;- 
axis anisotropy, arising from spin-orbit interaction in the 
2D QDs with energetically favorable heavy holes, 

O Nmu 

Hex = SziSzjSiVi -Rj), (10) 

i=l j=l 

where there are Nyin Mn spins at the position Rj. Here, 
Sz is the heavy-hole (pseudo)spin operator with projec¬ 
tions Sz = ±3/2, while Sz is the operator of the 2 - 
projection of the Mn-spin S = 5/2. Our theory does not 
include antiferromagnetic interactions between neighbor¬ 
ing Mn ions, which is relevant for QD’s doped with large 
Mn concentrations.®® 

Since Hex does not contain spin-flip processes, the to¬ 
tal wave function is a product of the hole and Mn-spin 
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parts.The partition function of the system can be cal¬ 
culated using a Gibbs canonical distribution, 

^ ( 11 ) 

which even for a single hole has a prohibitive complexity 
to be solved exactly. To calculate Z, in a typical QD 
with iVMn ~ 100 — 1000 and S=^I2 for Mn spin, one 
would need to solve replicas of the hole Schrddinger 
equation. 

To overcome this computational complexity and gain 
insight in magnetic ordering of open and closed-shell QDs 
we first use the MFA. We next consider a coarse-grained 
approach of discretizing the QD space to include spin 
fluctuations and examine the limitations of the MFA. We 
investigate the thermodynamics of the MP and the MBP 
formation and explore the finite-size effects by varying 
the number of magnetic impurity spins. We further cor¬ 
roborate our results and the influence of spin fluctuations 
using Monte Carlo simulations. 


III. MAGNETIC POLARONS (MPs) 

Early studies of MPs considered bulk magnetic semi¬ 
conductors where the localized carrier spin was provided 
by the donor or acceptor. por such bound magnetic 
polarons a finite extent of the carrier wave function leads 
to the alignment of only a small number of a nearby spins 
of magnetic impurities having many similarities with 
MPs in magnetic QDs. While problems of a conventional 
MFA in describing experiments on bound magnetic po¬ 
larons have been explained over thirty years ago,'^^ (spu¬ 
rious critical behavior and thermodynamic phase transi¬ 
tions in very small magnetic systems was removed after 
including spin fluctuations) such pitfalls continued to be 
repeated in describing magnetic QDs. 

We begin by considering a singly occupied QD [A/i = 1, 
Un), = 0, recall Eqs. (9) and (10)], the simplest real¬ 
ization of an open-shell QD,®^ and study the thermody¬ 
namics of the MP. We build the partition function by 
constructing a canonical Gibbs distribution, Eq. (11). 
The distribution function, Ds(iVfe,^fe), that describes the 
number of configurations of free spins in a given cell ex¬ 
pressed in terms of the microscopic parameter, (which 
can be viewed as a normalized magnetization), 

^s{.Nk. a) = E M ^ E ’ (12) 

{SjD V i=i ) 

where (recall Sec. 11) is the projection of the 

total Mn spin onto the z-axis, and the argument of the 
(5-function defines the normalization of ^ while the S- 
function is the distribution of Mn spins in a cell. We 
find fls{Nk,^k) oc exp[-Gs{^k)/kBT] (see Appendix A 
for details) with 

G 5 (Cfc) = kBTNk [^kBs^ i^k) - \nZs (Sg' (^ 0 )] , 

(13) 


being the free energy for non-interacting spins, where 
Zs{x) = sinh [(1 -I- l/2S')a::] /sinh [x/2S'], Bg^{y) is the 
inverse of the Brillouin function y = Bs{x).^‘^ 



FIG. 4. (Color online) (a) A magnetic QD doped with ran¬ 
dom paramagnetic Mn ions (green arrow) with one carrier 
spin density (red arrow), (b) The system lowers its energy 
through the exchange interaction and results in an antiferro¬ 
magnetic alignment between the hole spin and the Mn spins 
producing the MP. (c) The doubly degenerate QD energy level 
(1) splits with formation of the MP (2). The difference be¬ 
tween the nonmagnetic QD energy and the ground state en¬ 
ergy is the average exchange energy, Amp- 


A simple manifestation of a magnetic ordering in an 
open-shell QD occupied by a single carrier is the MP for¬ 
mation depicted in Fig. 4. Through exchange interactions 
between the carrier and Mn spins, once random param¬ 
agnetic Mn ions [Fig. 4(a)] acquire their spin alignment 
[Fig. 4(b)] (antiferromagetically coupled to a hole spin) 
and reduce the total energy of the carrier-Mn system. As 
a result of the hole-Mn exchange term in Eq. (10), the 
doubly degenerate heavy-hole energy level is split into 
two nondegenerate energy levels, as shown in Eig. 4(c). 
This corresponds to a red shift of the interband transi¬ 
tion energy as a function of time, which is observed in 
time-resolved photoluminescence experiments. 

In discrete space, for a given configuration of Mn spins, 
{^}, the two energy eigenvalues of Eq. (10) are, E± = 
±Amp/2, with spin-splitting energy 

23 

Amp = -;rB E ^kPMp{Rk)^k, (14) 

k 

where, 

PMp(i?fc) = 3l<).(i?fc)lV2, (15) 


is the heavy hole spin density at the fc’th cell, expressed 
in terms of the corresponding wave function cj). 

A MEA result for the spin-splitting energy A[[[p, is ob¬ 
tained by inter-relating the carrier and Mn spin densities 
(sz) and (^fe), respectively, in analogy of Eqs. (2) and (3), 




SkeT 


( 16 ) 
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and 


(&) 


Bs 


\ SksT 



(17) 


Substituting Eq. (16) and (17) into Eq. (14) gives 


aMF 

^MP 

2 


^Sj2NkPMpiRk)Bs f^^lff^tanh 


aMF ' 
^MP 

2kBT 

(18) 


In the unsaturated limit, A^p/2kBT -C 1, Eq. (18) gives 
a vanishing A^p at a critical temperature, 


r,T,7^MF 

'‘'B^C.MP 


P / 8(8 +1) 


R^kPwpiRk) 


1/2 


(19) 


The MP energy, Emp , is defined as the expectation value 
of Eq. (10). The MFA expression is 

aMF / aMF \ 


Having derived a ME theory description for the MP, 
one can make a connection between the MP and to that 
of a bulk nondegenerate magnetic semiconductors. Com¬ 
paring Eqs. (19) and (6), we see in both cases that the 
critical temperature has a linear dependence on (3. In a 
nondegenerate bulk semiconductor the concentration of 
donor/acceptor atoms is so small that the Pauli exclusion 
principle is ineffective and does not alter the carrier spin 
distribution which is described by the Boltzmann statis¬ 
tics. The magnetic impurity spins tend to align with one 
carrier spin and form a bound magnetic polaron.'^® The 
carrier spin in a nondegenerate magnetic semiconductor 
does not interact with other carriers spins. Thus, like the 
carrier in a magnetic QD, the carrier spins are free spins 
and will tend to align the spins of magnetic impurities. 

Next we turn our analysis to the MP energy. Mo¬ 
tivated by typical QD parameters,the blue/dashed 
line in Fig. 5 shows the ME behavior of the MP energy 
Eq. (20) as a function of T. For the chosen parameters, 
ME predicts a second order phase transition at charac¬ 
teristic temperature = 29 K. The phase transition 

occurs at a kBT comparable to the exchange splitting of 
hole levels. At this temperature both hole states of op¬ 
posite spin are approaching each other with equal proba¬ 
bility. Consequently, yielding vanishing average ME spin 
density and average exchange energy at a finite T. 

MFA neglects the possibility for the system to devi¬ 
ate from the minimum configuration, leading to 

phase transitions not allowed in nanoscale systems. To 
demonstrate the removal of the phase transition, we use 
the full partition function for the MP by summing over 
all configurations of The partition function is 


Zmp= E 

<T = ±1 \ B / 

( 21 ) 



FIG. 5. (Color online) A mean field solution for MP 
blue/dashed curve and the result obtained through the fluctu¬ 
ation approach red/solid. For simplicity, we approximate the 
hole wave function as uniform throughout the QD volume V, 
= 1/a/E.®® The parameters are NqP = —1.05 eV, the 
cation density, No = I/oq with ao — 6.1 A, E = Tvr'^hz, 
r = 5 nm, hz = 2.5 nm, Nmu ~ 90. 


with the spin index cr = ±1 for the heavy-hole spin Sz = 
|(T. Since the distribution Qs(Afe, ^k) is an even function 
of the integrals for cr = ±1 yield the same result and 

Zmb = exp ^ 2 fcBr^ ) ri ^s{Nk,^k)d^'^^. (22) 


Here the summation of cr = ±1 is done exactly, without 
neglecting the statistical correlation between cr and {^fc} 
as in the MFA. Now, using the steepest descent method 
as above, we have 


^MP 


’ c 




/ P8pMpiRk) 
V SfcsT 


Nk 


(23) 


The average exchange energy can be evaluated from 
Bmp = —kBTI3d{\n Zmp)/ dp as 


= -^8'^NkPMp{Rk)Bs 


where the corresponding results from Eq. (24) include 
the fluctuations of Mn spin and are compared in Fig. 5 
(red/solid) with the MFA results (blue/dashed). In many 
colloidal QDs the number of magnetic impurities is much 
smaller than used in Fig. 5 (Anhi = 90), enhancing the 
importance of the fluctuations and the corresponding dif¬ 
ference from the MFA solution. 

The MP shows different behaviors in the high- and 
low-T limits 

( -P for T ^ 0 

Emp oc < (25) 

[ -P'^/kBT for P8pMp{Rk)/kBT < 1, 

which correspond to saturated and unsaturated limits of 
magnetization, respectively. As shown in Fig. 5, Emp 
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has the l/T behavior for a wide range of temperatures. 
This can be understood as follows. As depicted in Fig. 6, 
a single carrier with uncompensated spin couples to a 
sum of many Mn spins. Therefore the carrier aligns with 
the majority of Mn spin, and a flip of an individual Mn 
spin would not affect the carrier spin, and consequently 
other Mn spins. For this reason, the Mn spins are in 
effect weakly interacting which results in a l/T Curie- 
like temperature dependence of Emp- 



FIG. 6. (Color online) Illustration of the finite-size effect for 
the MP (a) Spin density (blue/dashed) for fully aligned Mn 
spins (green). (b) A flipped Mn spin does not effect the carrier 
spin density (red/solid). 

We now show that the MP does not display a finite-size 
effect, when in a fixed QD volume we change the num¬ 
ber of Mn spins. Using the full MP partition function 
and the resulting Eq. (23), in Fig. 7(a) we show Tmp(T), 
normalized with respect to its fully saturated value for 
various number of Mn. Figure 7(b) shows the normal¬ 
ized Emp plotted at fixed T for various number of Mn 
ions. The normalized Emp remains constant. The finite 
size-effect was accurately predicted by the MFA. At the 
saturated and unsaturated limit, oc A^Mn- In Ap¬ 
pendix B we show how Emp will depend on a different 
choice of a carrier wave function. 

To summarize, we interpret the previous calculations 
in terms of the Ginzburg-Landau approach of phase tran¬ 
sitions. While the approach is not reliable due to the 
nanoscale size of our system, we may gain important 
insights. Typically the partition function is described 
by a functional integral over the magnetization ^ [see 
Eq. (21)]. In the low-T limit (or |AMp|/2fcBT ^ I) with 
gAMp/s/csr _j_g—AMp/zfcBr pg g|AMp|/ 2 fcBT gq (21) one 

can write the free energy 

FMp{0^-gi\^\-kBTS{0. (26) 

with the entropy S'(^) even in Tmp(0 is depicted in 
Fig. 3(a). However, the finite potential barrier separat¬ 
ing the degenerate minima at ±^min does not prevent the 
thermal fluctuations between the local minima. There¬ 
fore the correct solution (^) = 0 is not predicted by the 
MFA to Tmp(C). 

We consider the Ginzburg-Landau approach defined 
with a different variable, the observable quantity of the 
exchange energy X = E^^ = crA mp/ 2 [see Eq. (14)] with 
the carrier spin index cr = ±1. The linear dependence in 


^ originates from the hnite carrier spin of the open shell. 
Since the carrier spin does not contribute to the entropy, 
the entropy in X is directly related to S{^), and we derive 
the free energy in the order parameter X = Tgx as 

EMPiX)^X-TSiX), (27) 

as depicted in Fig. 3(b). Again, the entropy S{X) is 
an even function of X and Eq. (I) results for the mag¬ 
netic polaron. Unlike Tmp(C), T'mp(A) possesses only 
one global minimum at a negative hnite X and {X) < 0 
at all T, and the mean-held interpretation of Eq. (27) 
gives a qualitatively correct prediction. The thermody¬ 
namic solution of a hnite (A) manifests itself in the hnite- 
scale independence in {X)/NMn in Fig. 1(d). 



FIG. 7. (Color online) (a) The MP average exchange energy 
normalized by its T = 0 K value for 10 Mn (red/solid) and 
20 Mn (black/dashed) as function of T. (b) Normalized MP 
average exchange energy as a function of the total number of 
Mn at fixed T. 


IV. MAGNETIC BIPOLARONS (MBPs) 

We next turn to a magnetic QD containing two holes 
[V/^ = 2, U 2 = U, Eqs. (9) and (10)]. Closed-shell 
fermionic systems, such as noble gases, are known for 
their stability and the total spin-zero ground state, mak¬ 
ing them magnetically inert. Thus it would seem that 
this simple example of a two-hole closed-shell QD doped 
with Mn would not allow magnetic ordering. However, 
the Mn doping does alter the magnetic properties of 
closed-shell QDs. The corresponding ground state, which 
is neither a singlet nor a triplet, allows ordering of Mn 
spins, owing to the spontaneously broken time-reversal 
symmetry.®® 

To lower the hole-Mn system energy through exchange 
interaction, there needs to be a nonvanishing hole spin 
density. During the MBP formation, an initially random 
Mn spin orientation [Fig. 8(a)], in the presence of two 
holes acquires a spin alignment [Fig. 8(b)]. The emer¬ 
gence of a nonvanishing local hole spin density, while the 
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total hole spin density remains zero, is characteristic for 
a spin pseudosinglet,®’®® which can be understood from 
a simple perturbation picture. The exchange interaction 
admixes higher (single particle) orbitals to the ground- 
state s orbital. Specifically, as shown in the Fig. 8(c), 
the mixing of s and px orbitals leads to the spin-Wigner 
molecule,® a spin-analog of the Wigner molecule.®^“®® In 
contrast to the Wigner molecule, where the spatial carrier 
separation originates from the Coulomb repulsion,®® here 
the dominant contribution of such separation is typically 
the carrier-Mn exchange energy. 



Magnetic Bipolaron Formation 


(C) 


Mixing s and p orbital 


b:) 

O 




Ep 

Es 



FIG. 8. (Color online) (a) A magnetic QD doped with ran¬ 
dom paramagnetic Mn ions spins (green arrows). (b) The sys¬ 
tem lowers its energy through the exchange interaction lead¬ 
ing to a nonvanishing spin density and formation of the MBP. 
Red arrows and lines in (a) and (b) show how the two-carrier 
spin density changes due to the presence Mn. (c) The doubly 
occupied QD energy level (1) lowers its energy through the 
formation of the MBP (2). The difference between the non¬ 
magnetic QD energy and the MBP ground state energy is the 
average exchange energy Sex- 


taken to be constant, and the third term is the average 
exchange energy between the MBP spin density, 

6 ^ 

PMBpiRj,£x) = -—^s{Rj)px{Rj), (30) 

at the site of Mn. For a nonmagnetic system, the two- 
particle energy for the ground state is Eg = 2hbj and for 
the p-state is Ep = ihuj. To get close to the eigenstate of 
the system, we seek the that minimizes Eq. (29) 


_ Ambp 

Epsji -h + A^gp 


(31) 


and minimized energy, 

Emin = 2 (Es + Ep — ■\JEpg 4A^gp^ -h C/, (32) 

where Epg = Ep — Eg is the energy difference between 
the s and Px orbital, and the spin splitting due to orbital 
polarization. 


Ambp 


AfMn Nc 

/3 ^ s{R,)px{Rj)S,,x = Epg ^ xk^k, (33) 

j=l k=l 


with a dimensionless magnetic susceptibility like term 


Xk 


SNkPs{Rk)px{Rk) 

E 

-E-^pS 


(34) 


To obtain the full partition function we need to sum 
over all configurations of ^k- Since we are in the fepE <g; 
Epg regime, we consider only the spin “singlet” ground 
state and neglect the possibility of a triplet state.®® The 
partition function is 


The corresponding pseudosinglet wave function®® at 
each snapshot of Mn conhguration, {^k}, is 

N 

(t>ps{ri,r2) = [u{ri)d{r2) |t;4-) - u{r2)d{ri) |4.;t)] , 

(28) 

where N = 1/(1 -|- e®) is the normalization constant, 
u{r,{^k}) = s{r) + £x{{^k})Px{r), d{r,{^k}) = s{r) - 
£xi{^k})Px{r), and ex({Cfc}) is a mixing (variational) pa¬ 
rameter that depends on In choosing our varia¬ 

tional wave function we neglect overlaps between like- 
orbitals, (s — s and px — Px), which results in the loss of 
fluctuations in the total magnetization at the site of Mn. 
These fluctuations are small in the MBP regime. 

The variational energy E{ex) is 

E = f E PMBp{R,.ex)S,,x. (29) 

E “T Srp O . 

^ J =1 

where the first term is the sum of the kinetic and po¬ 
tential energy, the second term is the Coulomb energy, 


Embp = / Y[ ns{Nk,^k)d^^^. (35) 

For simplicity, to evaluate Eq. (35), we use a two site 
model by placing Nk Mn at two opposite sites equally 
spaced from the origin along the x-axis. The sites are 
chosen such that s{Ri)px{Ri) = —s{R2)px{R2), where 
El and R 2 are the position of the Mn ions at site 1 and 
site 2, respectively. For this two site problem, Eq. (33) 
reduces to Ambp = EpsX^_, where 

(36) 

and y = |yi| = 1x21, Ni= N 2 = AMn/2. 

To investigate the MBP as a function of E, we begin by 
using the MFA. It can be shown that the minimum of the 
free energy, Embp(? -) — Emin (a,6)+Gs(a)+Gs(6), 
lies on the = —^2 line. The free energy becomes 

Embp(^-) = Emi„(C-) + 2Gs(e-/2). (37) 














Minimization of the MBP free energy, Eq. (37), with re¬ 
spect to gives a self-consistent equation for A|^gp, 


A MF 
^MBP 


— ^psX^S 


^xAmbp 


NMnkBT^El^+AiA^l^r) ’ 

(38) 

where Ajjjgp = EpsX^^^■ In the unsaturated limit, 
^MBp/^ bT < 1, Eq. (38) gives a vanishing at 

a critical temperature, 




8S{S+l)x^Ep 

^Nun 


cx 


(39) 


Unlike the MP, Eq. (19), the MF critical temperature 
for MBP is quadratic in /?, through x^, a result similar 
to that of degenerate holes in bulk DMS, Eq. (8). 




to give 


E. 


MF 



(40) 


The green/dashed line in Fig. 10(a) shows the T- 
dependence of Eq. (40), while Fig. 10(b) shows the mean 
field behavior of the averaged product of the normal¬ 
ized magnetization at the two sites, (m 2 mi). There is 
an antiferromagnetic correlation between the product of 
normalized magnetization at the two sites. As a conse¬ 
quence of thermal spin fluctuations, the average exchange 
energy and the magnitude of (m 2 mi) decreases and with 
increasing T eventually vanishing when the MFA for 
the specific parameter set yields a vanishing carrier spin 
density resulting in a second order phase transition at 
T^mf _ 1 c 
-‘C.MBP — 

Exact integration of the MBP partition function 
Eq. (35), which correctly includes spin fluctuations, is 
needed to remove the phase transitions predicted by MF 
theory. The average exchange energy is obtained through 


Afix — 


1 


^MBP _ 

9Emin(Cli C 2 ) 


13- 


dp 


d^i I d^2^siNi,^,)nsiN2,^2) 

-F'min 


A T-dependence of this average exchange energy is 


FIG. 9. (Color online) Illustration of the finite-size effect for 
the MBP (a) Spin density (blue/dashed) for fully aligned Mn 
spins (green), (b) A flipped Mn spin changes the carriers spin 
density (red/solid). 

We see that pMBp(ea;) in Eq. (30) is coupled linearly to 
the magnetic ordering of the Mn spins via Eqs. (31) and 
(33), a result that is analogous to Pauli paramagnetism 
where the magnetization of a free electron gas is propor¬ 
tional to the strength of an external magnetic field. Thus, 
any small change in the Mn configuration will result in 
a linear response in the distribution of the carriers’ spin 
density. The strength of the response is determined by a 
Pauli-like susceptibly term given by Eq. (34). Figure 9(a) 
shows that one Mn spin “sees” the carrier orbital spin 
while one carrier spin “sees” all Mn collectively. Since 
the height of the spin density is determined by , which 
is dependent on the configuration of Mn spins, if one Mn 
spin changes. Fig. 9(b), the exchange field arising from 
the Mn field is strong that the carrier responds to the 
change. Therefore, inverting one Mn spin, Fig. 9(b), will 
change the amplitude of the spin density. Through the 
exchange interaction between the holes’ spin density and 
Mn spins, all other Mn will respond to the changing spin 
density. Thus, the Mn are indirectly interacting with one 
another through the hole. 

The average exchange energy for the MBP is obtained 
the same way as the MP, using the solutions of Eq. (38) 



FIG. 10. (Color online) (a) T-dependence of the average 
exchange energy and the removal of the phase transition 
for the two site model of the MBP, [Eq. (40)] with a non¬ 
fluctuating carrier spin density (green/dashed) and a fluctu¬ 
ating spin density (red/solid) [Eq. (41)]. (b) T-dependence 
of the product of magnetization of each site, (m 2 mi), for a 
fluctuating spin density (red/solid) and nonfluctuating spin 
density (green/dashed). The parameters are hz — 2.5 nm, 
tko = 30 meV, ml — 0.21, NqP — —1.05 eV, U = 30 meV, 
i?i ,2 = ±2 nm and Nmu = 10. 

shown by the red/solid lines in Fig. 10(a). With the 
inclusion of spin fluctuations, the phase transition is re¬ 
moved resulting in a finite Egx at T > Tc- The product 
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of the normalized magnetization of each site was calcu¬ 
lated numerically to demonstrate the antiferromagnetic 
correlation between the Mn spins at the two sites. 

Monte Carlo simulations (see appendix C for details) 
verify the theoretical prediction for the MBP. s — p lev¬ 
els coupled via Mn spins are directly diagonalized and 
the variational assumption, Eq. (28), was not used. Mn 
ions (5 per site) were placed 2 nm from the center of 
the QD along the i-axis. We include s, Px and Py or¬ 
bitals in the simulation. The Coulomb energy, as before, 
is a constant. The blue dots in Fig. 10(a) shows the T- 
dependence of the average exchange energy and the prod¬ 
uct of the magnetization per site is shown in Fig. 10(b). 
There is an excellent agreement between the fluctuation 
approach and Monte Carlo simulations. Furthermore, 
the antiferromagnetic correlation between the two sites 
demonstrates the formation of the MBP. 

As we did for the MP, we investigate two limiting cases 
for Aex. In the saturated regime, the Mn are maximally 
aligned at each site with their spin pointing in opposite 
directions with |^i — ^2! —t 2. In this limit, the MF 
expression, Eq. (40), is valid. In the opposite unsatu¬ 
rated high-T limit, fa 1 can be taken in 

Eq. (41). Finally, we obtain the limiting behavior of Fex 
within the regime x -C 1, 

-G^NunT^J^BP « (^Mn^)" for T ^ 0 

MBP oc for T -)> oo 

(42) 

This shows that, unlike the MP, the coupling dependence 
is always in /3^. This is due to the fact that the carrier 
spin density is polarized from zero by the coupling to Mn 
spins, unlike the uncompensated carrier spin in the MP. 

The MBP displays a finite-size effect different from 
the MP. We plot the T-dependence of the average ex¬ 
change energy, Eq. (41), normalized by its low-T value 
in Fig. 11(a) with a varying number of Mn in a given 
cell. As the number of Mn increases, the normalized Fex 
fluctuation tail decays toward the MF solution. This is 
demonstrated in Fig. 11(b) at a fixed ratio of T/Nmu- 
As the number of Mn ion spins increases toward the 
thermodynamic limit, the normalized exchange energy 
decays. This size dependence is expected from the MF 
result where a phase transition to zero exchange energy 
occurs. The l/A^Mn dependence in Fex reflects the ther¬ 
modynamic limit towards the vanishing order parameter. 

To summarize, like in the MP case, the MFA can be 
used to understand the statistical properties of the MBP. 
Applying the variational treatment to the carrier spins, 
we express the exchange energy as Aex = —a(Ci — ^2)^ 
with the Mn spins ^1^2 at position 1 and 2. This quadratic 
dependence (.^1 — ^2)^ originates from the carrier spin 
being linearly polarized out of the spin-singlet in closed- 
shell systems. With the variable X = — ^2, the free 

energy becomes 

(43) 



FIG. 11. (Color online) The finite-size effect in a doubly oc¬ 
cupied QD. (a) T-dependence of the normalized MBP average 
exchange energy for 10 Mn (red/solid), 20 Mn (black/dashed) 
and 40 Mn (blue/dash-dotted), thus 5, 10, and 20 Mn per site, 
respectively. Mean-field solution with 10 Mn (black/dotted) 
and 40 Mn (green/dashed), (b) The normalized mean-field 
solution with varying number of Mn for a fixed T/Nmu- 


Confining the discussions to unpolarized states (.^1 + 
^2 = 0), the entropy becomes an even function of X and 
Eq. (1) results with 51 = 0 [see Fig. 1(b)]. We note that, 
in contrast to the Eq. 1, here the first X'^ term in the 
free energy is not of an entropic origin. In a finite sys¬ 
tem {X) = 0 and the exchange energy Fex = 
remains finite. While the MFA incorrectly predicts a 
phase transition to = 0 beyond a phase transition 
temperature, it correctly justifies the finite-size scaling 
limit Eex/Nun —f 0 as A^Mn —>■ 00 as shown in Fig. 1(d), 
in a sharp contrast to the MP case. The open or closed- 
shell electronic structures leads to fundamentally differ¬ 
ent statistical properties in the quantum dot magnetism. 


V. CONCLUSIONS 

Many of our findings for magnetic polarons and bipo- 
larons can be also applied to higher carrier occupancy 
of quantum dots where it is important to distinguish 
if they form open- or closed-shell systems which leads 
to a qualitatively different classes of magnetic ordering. 
Performing a thermodynamic analysis in these nanoscale 
magnets, we reveal the limitations of a mean-field ap¬ 
proximation and the necessity for a more accurate theo¬ 
retical framework that would correctly include spin fluc¬ 
tuations. Our results show that a careful choice of the 
order parameter in the mean-field approximation (using 
the exchange energy, rather than magnetization) removes 
spurious phase transitions for magnetic polarons, but not 
for magnetic bipolarons. In the later case, the phase tran¬ 
sitions are removed by including spin fluctuations within 
the coarse-grained method and Monte Carlo simulations. 

The conventional mean-field theory, known from bulk 


A'mbp(^i,C2) ~ —aX^ — 
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systems with nondegenerate or degenerate carrier den¬ 
sity, reveals important differences between the magiretic 
polarons and bipolarons. Surprisingly, we can introduce 
a very simple meair-field form of a free energy to ac¬ 
curately describe qualitatively different finite-size effects 
and distinct thermodynamic limits in magnetic polarons 
and bipolarons with the change of the number of mag¬ 
netic impurity spins. These findings remain unchanged 
once we carefully include spin fluctuations, further jus¬ 
tifying our simple description and a pictorial difference 
between the finite-size effects in magnetic polarons and 
bipolarons, Figs. 6 and 9, respectively. 

Similar to our prediction for an unexpected thermally- 
enhanced magnetic ordering in quantum dots,^^ a judi¬ 
cious use of the mean-field approximation and awareness 
of its artifacts could provide important insights in unex¬ 
plored aspects of nanoscale magnets. For example, we 
expect that the mean-field description of the different 
finite-size scaling in magnetic order discussed for mag¬ 
netic polaron and bipolaron will also apply to other open- 
and closed shell quantum dots with higher carrier occu¬ 
pancy. A mean-field calculation of the critical tempera¬ 
ture could also reveal a different power-law dependence in 
the exchange coupling constant for the exchange energy 
of open- and closed-shell systems. 

In contrast to magnetic polarons, much less is known 
about magnetism in closed-shell systems, often simply 
implying that the magnetic ordering is completely ab¬ 
sent. Therefore, to test our predictions for magnetic 
bipolarons it would be important to focus on the ex¬ 
perimental realization of multiple carrier occupancy in 
quantum dots. The simple creation of excitons is not 
sufficient. A simultaneous presence of single electron 
and hole effectively just renormalizes the exchange cou¬ 
pling with magnetic impurities of magnetic polarons. 
Instead, photoexcitiation, using chemical and electro¬ 
static doping,®^ should create a pair of holes or electrons. 

Another possibility would be to fabricate quantum 
dots from novel Mn-doped II-II-V dilute magnetic semi¬ 
conductors. These systems provide an independent 
charge and spin doping and would therefore be suitable to 
test formation of nanoscale magnetism for a wide range 
of parameters.They share with (II,Mn)VIs an iso- 
valent character of Mn-doping, removing the solubility 
constraint of (III,Mn)Vs an obstacle for fabricating mag¬ 
netic quantum dots.^° Unlike (II,Mn)VI, Mn-doped II-II- 
V systems can separately attain different carrier densities 
through independent charge doping and thus readily al¬ 
ter the strength of the Mn-carrier exchange coupling. 
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Appendix A: Mn Distribution Function 


We discretize the QD space with magnetic mo¬ 
ments in a given cell. The distribution function is 
found by summing over all configurations of the Mn 
in a given cell and is determined by Eq. (12). Us¬ 
ing the integral representation of the 5-function, 5(a;) = 
(l/27r) the partition function is 


ns(N,A) = f 




(Al) 

Introducing a complex variable h = iXS^ Eq. (Al) be¬ 
comes 


ns{Nk,^k) = ^ n 

27rz J —200 


(A2) 


where Zs{x) = sinh [(1-|-l/25')2:] / sinh [a;/25']. The in¬ 
tegrand in Eq. (A2 ) is sharply peaked (Gaussian-like), 
therefore we can approximate Eq. (A2) by performing 
the method of steepest descent. We deform the contour 
in the complex plane to pass through a saddle point in 
the direction of steepest descent. By Taylor expansion of 
the function in the exponent of Eq. (A2) and performing 
a Gaussian integral over h we obtain 

where x(h) = d'^\nZs{h)/dh‘^ and Gs{^k,T) is the 
Gibbs free energy, recall Eq. (13), obtained through a 
Legendre transformation. 


Appendix B: Magnetic Polaron Renormalization 


Since the MP is not influenced by the finite-size effect 
(recall Figs. 6 and 7), we can employ a variational ap¬ 
proach to study the influence of the wave function renor¬ 
malization on the MP properties in a QD. Starting from 
the MP partition function in Eq. (23), we approximate 
the MP wave function by a single s-orbital of a 2D har¬ 
monic oscillator, constant over the the QD height, 


4>{r) = l/{^/h^LMp)e (b1 ) 


Using a variational approach^^’^^ we determine the 
width, Amp, that minimizes the total free energy func¬ 
tional, Fmp- The MP free energy functional is given by 


Fmp = 


hu 


t2 
Nmti 

- ksT In 


/■2 

Ll 


- kBTln2 


Zs 


f ldSpMp{Rj) 

V 3fc^T 


(B2) 
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where the first two terms are the sum of the kinetic and 
potential energy with Lq = the third term is 

from the hole spin degeneracy, and the final term is due to 
the exchange interaction of the hole spin density pMp{r) 
at the site of Mn spins, Rj. As an approximation, we 
consider a homogeneous distribution of Mn, and trans¬ 
form Nqxmh f (PR (the continuous limit) where 

XMn is the Mn fraction per cation and Nq is the density 
of cation sites. The average exchange interaction for the 
MP, Amp, can be derived to obtain 

Amp = J d^RpMp{R)Bs > (B3) 

where A^ax = a^MnlAo/3|5'. Amp is found by numeri¬ 
cally minimizing Eq. (B2) to obtain the most probable 
width, Amp- This width determines pmp by combining 
Eqs. (Bl) and (15) and yields Amp from Eq. (B3). 

It is instructive to now compare how various forms 
of the carrier wave function affect the T-dependence of 
the Amp, shown in Fig. 12(a). We choose XMn = 2.6%, 
hz = 2.5 nm, fiw = 30 meV, m)) = 0.21, Lq = 3.5 nm- 
the characteristic width in the absence of Mn spins, and 
Aq/? = —1.05 eV, which within the classical radius of 
this harmonic confinement yields Amh = 90, as in Fig. 5. 
From Fig. 12(a), where we compare our results for the 
variationally obtained Amp(T) with the wave function of 
a fixed width at Amp(T) = Lq, we see that the wave 
function renormalization have a very small influence on 
Amp (A). In fact, both of them are very similar to the 
variational Amp (A) for a constant wave function in Fig. 5. 



FIG. 12. (Color online) (a) T-dependence of the MP en¬ 
ergy [Eq. (B3)] with the variationally-obtained (“renormal¬ 
ized”) (red/solid) and the fixed wave function for a harmonic 
QD confinement without Mn spins, (b) The corresponding 
T-dependence of the most probable width for a renormal¬ 
ized wave function, Tmp (red/solid) and fixed wave function 
(blue/dashed), Tmp = To, from Eq. (Bl). (c) T-dependence 
of the Tmp from a variational method (red/solid) and Monte 
Carlo simulations (blue/dotted). 

From Fig. 12(b) we see that the wave function renor¬ 
malization itself is a small effect. The red/solid curve 
shows how the wave function width, normalized to the 


nonmagnetic width, Lq, varies with A. As the A in¬ 
creases from A = 0 K, the Mn spins coupled to the tail 
of the wave function are more prone to thermal excita¬ 
tion. The wave function shrinks to attain a more energet¬ 
ically favorable configuration by increasing the exchange 
energy gain from the polarized Mn spins near the cen¬ 
ter of the QD. Eventually, thermal excitation overcomes 
the magnetic energy, and the system relaxes continuously 
to a nonmagnetic state resulting in Tmp = Lq at large 
A, as there is no energy gain from the wave function 
renormalization. From the variational approach we see 
an additional MP localization: Tmp < Lq, while the non¬ 
monotonic Tmp (A) implies also a nonmonotonic effective 
exchange field^^ due to the MP formation. 

To further verify the renormalization effect of the MP 
wave function, we implement Monte Carlo simulations, 
removing the need of a variational calculations of the 
wave function. To approximate the MP wave function we 
include s, p, and d orbitals in our simulation and allow for 
mixing of these orbitals. In Fig. 12(c) we see that both 
variational (red/solid) and Monte Carlo (blue/dotted) re¬ 
sults agree well with each other and shown again a non¬ 
monotonic Tmp(A), noted also in Fig. 12(b). The slightly 
smaller wave function renormalization in Fig. 12(c), as 
compared to that in Fig. 12(b), is a consequence of fewer 
Mn spins in the middle region of the QD. 

Appendix C: Monte Carlo Simulations 

Monte Carlo simulations were used to approximate so¬ 
lutions to the Schrodinger equation 

A({A})|$)=A({A})|$), (Cl) 

for a fixed finite orthonormal basis |$) at a given Mn spin 
({S'}) configuration. The calculation entails guessing a 
Mn configuration at a given A, producing a matrix rep¬ 
resentation of H{{S}) in a finite basis, and solving the 
eigenvalue problem. 

The calculation begins by defining a 2D harmonic QD 
and solve for single heavy hole QD levels and eigenfunc¬ 
tions without any Mn atoms. We truncate the states up 
to the first A-orbitals. We obtain A noninteracting wave 
functions energy A° with a For a given 

configuration of Amh Mn spins, {S^i, Sz2, •••, SzNMn}^ ^e 
construct a (2A) x (2A) matrix, 

H = J2E° \4>L) (<(>11 + E 9nn'a{{Sz}) (<(>11 ■ 

ncr nn' 

(C2) 

The interaction constant g is 

yO ^Mn 

gnn'a{{Sz}) = f (s.,.S., ) <)>1 (R, ), 

(C3) 

where Sz is the spin of the carrier and Rj is the posi¬ 
tion of the Mn ion. We diagonalize the single heavy hole 
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Hamiltonian at a snapshot of {5’^} and obtain eigenvalues onalize H{{S'^}) and obtain new eigenvalues iS'dS''}). 
Ei{{Sz}) and eigenvectors Ci^na{{Sz})- We then propose 
a different Mn-conhguration, {5'd:diag- 
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